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Abstract. We consider two-dimensional (2D) arrays of self-organized semiconductor 
quantum dots (QDs) strongly interacting with electromagnetic field in the regime of 
Rabi oscillations. The QD array built of two-level states is modelled by two coupled 
systems of discrete nonlinear Schrodinger equations. Localized modes in the form 
of single-peaked fundamental and vortical stationary Rabi solitons and self-trapped 
breathers have been found. The results for the stability, mobility and radiative 
properties of the Rabi modes suggest a concept of a self-assembled 2D soliton-based 
nano-antenna, which should be stable against imperfections In particular, we discuss 
the implementation of such a nano-antenna in the form of surface plasmon solitons 
in graphene, and illustrate possibilities to control their operation by means of optical 
tools. 
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1. Introduction 

Experimental and theoretical studies of solitons in nonlinear lattices are the field of 
intensive activity in several branches of physics. Its large part is focused on nonlinear 
optics [1] and matter waves (i.e., Bose-Einstein condensates, BEC). The latter may be 
treated in the mean-field approximation [2], or in terms of the quantum Bose-Hubbard 
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model [3, 4] (see book [5] for a systematic summary). In these contexts, the discrete 
nonlinear Schrodinger (DNLS) equation, and systems of the coupled DNLS equations, 
are ubiquitous dynamical models for the description of nonlinear lattices [6]. In the 
framework of these studies, the existence, stability, and dynamics (including mobility) 
of discrete self-trapped modes (discrete solitons) is a topic which has drawn a great deal 
of interest. 

In multi-component systems, couplings between components may be linear, 
nonlinear, or both linear and nonlinear [6]. In the context of optics, systems of 
linearly coupled DNLS equations are relevant to various applications. In particular, 
the linear coupling between two polarization modes in each core of a waveguide array 
may be induced by a twist of the core (for linear polarizations), or by the birefringence 
(for circular polarizations). On the other hand, in the BEC the linear coupling may 
be imposed by an external microwave or radio-frequency field, which can drive Rabi 
oscillations [7, 8] or Josephson oscillations [9, 10] between two boson populations. In 
terms of the dynamical analysis, an essential issue is the possibility of the spontaneous 
symmetry breaking in linearly coupled discrete systems [11], 

The progress in the fabrication of nano-scale electric circuits [12], lasers [13], 
waveguides and antennas, for microwave, tera-Hertz and optical frequency ranges, is 
another strong incentive stimulating studies of nonlinear lattices. Particularly significant 
nano-dimensional elements are self-organized quantum-dots (QDs) [14] embedded into 
semiconductor dielectric host media. The possibility of strong interactions of QDs with 
optical fields gives rise to Rabi oscillations (RO) between electron-hole populations of 
ground and excited states. The basic RO model amounts to a two-level atom strongly 
coupled to an external electromagnetic field (the Jaynes - Cummings model) [15]. For 
more complex structures, such as large molecules [16] or QDs, this model has been 
modified to account for additional factors: anisotropy [16], local-field corrections [17]- 
[20], broken inversion symmetry [21], etc. In these settings, crucially significant are 
collective coherent effects, stipulated by the inter-dot coupling inside the QD array. 
Oscillations of the QD population between levels in an isolated two-level system may 
be considered as energy exchange between the system and the ambient electromagnetic 
field. On the other hand, particle transfer (e.g., inter-dot electron, or electron-hole 
tunnelling) leads to the exchange of the quasi-momentum between charge carriers and 
photons. These mechanisms govern the spatiotemporal RO dynamics in the QD chains, 
in the form of the propagation of travelling RO waves and wave packets (Rabi waves) 
[22]-[24], 

In the QDs arrays the transport can also be realized by tunnelling via the inter¬ 
dot dipole-dipole coupling, such as Forster interactions [14], and by the radiation-field 
transfer. In Ref. [23], the tunnelling was assumed to be a dominant mechanism of the 
inter-dot coupling. In that case, a significant role is played by local-field effects in the 
QDs, even in the weak-coupling regime [17, 18]. The local-field effects, enhanced by the 
strong light-QD coupling, modify the conventional RO dynamics [22]-[24] . As a result, 
specific nonlinear terms appear in equations of electron-hole oscillations inside the QDs, 
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braking the superposition principle for the single-particle wave function. The role of 
local fields in the formation of excitonic RO in self-assembling QDs was experimentally 
studied in Refs. [25, 26]. 

The nonlinearity affects the Rabi-wave propagation too. As was shown in Ref. 
[27], solitons and self-trapped breathers in the one-dimensional (ID) QD chains can be 
created with the help of the nonlinearity. The model developed in Ref. [27] explores the 
dynamics of the probability amplitudes of the ground and excited states of QDs, with 
included nonlinear cross-phase-modulation terms and linear coupling between the two 
states. The phase shift between constants of the inter-site hopping is a new property 
introduced geometrically in that ID setting. 

Recently, theoretical and experimental studies of nano-antennas have been in a 
focus of work performed by many groups. The antenna is defined as a device which 
transforms the near field into the far field, or vice versa (transmitting and receiving 
antennas, respectively). Transmitting devices produce the spherical wave in a remote 
spatial area, with electric field 

,2 0 — iuj{t—R/c ) 


ERad — 


e e or e 

-n, 0 - 


-F(9, tp) + c.c., 


( 1 ) 
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where R , 0. ip are the spherical coordinates, with the origin set at the center of 
the antenna, F(9, <p) is a normalized angular radiation pattern, e 0 is the vacuum 
permittivity, c is the free-space light velocity, eg is the unit vector along the local 
direction of coordinate 9, with E 0 and uj being the amplitude and frequency of the 
electric field, respectively. The purport of using the antennas is their ability to provide an 
interface between local information processing, which uses electrical signals, and the free- 
space wireless transmission of data encoded in various parameters of the electromagnetic 
waves, such as the amplitude, phase and frequency. 

Directional antenna properties (the dependence of the emitted field on both the 
azimuthal and polar angles) are characterized by the radiation pattern. It is determined 
by the near-field distribution produced by the signal source placed inside the antenna. 
It depends on the frequency, antenna configuration and geometric parameters. An 
important species of the device is represented by phased antenna arrays, i.e., systems of 
a large number of identical emitters with a phase shift between adjacent ones [28]. The 
progress in the fabrication of nano-antennas, see Refs. [29]-[31] for review, manifests the 
general trend of implementing radio-communication principles in the optical-frequency 
range. In particular, it stimulates bridging the realms of macroscopic antennas and nano¬ 
antennas, using new materials with specific electronic properties. In this context, it is 
relevant to mention carbon nanotubes and nanotube arrays [32]-[35], plasmonic noble- 
metal wires [36, 37], graphene nanoribbons [38], semiconductor QDs [39, 41], etc. These 
types of nano-antennas have been offered as promising elements for industrial design and 
commercial manufacturing. However, many limitations for their use are stipulated by 
difficulties of their operational control - for instance, rotation of the radiation pattern 
by an electric or optical drive. This difficulty, along with the growing interest to the 
field, motivate the search of new physical principles for fabrication and tuning of nano- 
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antennas. A promising possibility for that is suggested by using strong nonlinearities in 
nanostructures (in particular, remarkable nonlinear properties of graphene [42, 43]). 

In this work, we propose a previously unexplored principle for the realization of 
nano-antenna arrays in the form of self-trapping discrete solitons in 2D nonlinear lattices 
built of semiconductor QDs. In this context, dynamically stable soliton structures 
(fundamental discrete solitons, vortex solitons, and breathers) may provide promising 
mechanisms of controllable creation of the antennas with various geometric shapes. 
The model for the 2D array of QDs can be developed from its ID counterpart which 
was elaborated in Ref. [27]; obviously, the transition from ID to the 2D geometry 
is a crucially important step for modelling antennas which emit in the plane of the 
semiconductor structure. We demonstrate that different types of the global drive, 
viz., plane and cylindrical waves, can excite Rabi oscillations and waves in the system. 
Using known numerical techniques developed for the study of discrete solitons [6], we 
construct basic localized modes of the corresponding 2D DNLS equations and test their 
stability and mobility. In addition, we demonstrate how discrete Rabi solitons induce the 
dielectric polarization, and thus build the near-field of the soliton-based nano-antenna. 
On the contrary to other radiation zones, the near field keeps itself in a quantum state, 
which helps to establish the correlation of quantum effects and field transformation 
between different spatial zones. We also present the far-field structure produced by the 
given near field, which takes into account the quantum nature of the underlying Rabi 
oscillator. 

The paper is organized as follows. In Sec. II, we introduce the model for the 2D 
QD array interacting with the driving electromagnetic field. As mentioned above, the 
plane- and cylindrical-wave drives are considered. The nonlinearity is induced as the 
local-field effect. The model is analyzed numerically, with the purpose to reveal the RO 
dynamics in the QD arrays. In Sec. Ill, necessary technical ingredients are introduced, 
including the dispersion relation for the linearized system, and general expressions for 
radiation properties generated by solitons in the 2D QD array. In Sec. IV, results are 
presented in detail, and the relation of the present model to 2D nano-antenna arrays 
is proposed. The important issue of maintaining the coherence of emitters building 
the soliton nano-antenna is considered in Sec. IV too, and possibilities for operational 
control of the device by means of optical tools are outlined. The paper is summarized 
in Sec. V. 

2. Two-dimensional model equations 

2.1. The plane-wave drive 

We consider finite 2D rectangular N\ x N 2 array of identical QDs with the square unit 
cell of size a x a (a is the inter-dot distance). The position of QDs is defined by a pair 
of discrete indices, p = —N\/ 2,..., N\/2 and q = —N 2 / 2,..., N 2 / 2. The array is exposed 
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to the classical travelling optical wave of in the xy- plane with electric field 

E(x, y, t)=Re {E 0 exp [ik (a: cos do + y sin do) — loot]} , (2) 

where k = k(u>) is the wavenumber and d 0 is the propagation angle (see Fig. 1 (a)). The 
driving field may represent different types of the electromagnetic modes, the first one 
being a plane wave with a nonzero E z component. Another example is a surface plasmon 
guided by the plane boundary between the noble metal and transparent medium [14], on 
which the QD array is installed. The latter example comes with the frequency dispersion 
of a plasma medium, which manifests itself in a nonlinear dependence k = k(cu). We 
take QDs as identical two-level dissipationless systems with the transition energy, hu>o, 
corresponding to the transition between excited and ground-state electron orbitals, \a P:q ) 
and | b Pt q). The transition dipolar moment is p = e z fi, and only intraband transitions are 
taken into account [23]. Each QD is coupled by the electron tunnelling to four nearest 
neighbors in the 2D array. 



Figure 1 . Schematically presented plane wave (a) and cylindrical wave (b) drive of 
2D QD array placed in the (x, y ) plane. Dotted curves display wave fronts in both 
cases, k is the wave-vector, 9q is the propagation angle and a is the inter-dot distance. 


It is adopted that the light interaction within the QD chain takes place in the 
resonant regime, i.e., the frequency detuning is small in comparison to both the light 
and quantum-transition frequencies. Following the rotation-wave-approximation [15], 
we omit rapidly oscillating terms in the equations of motion. These natural assumptions 
follows the Rabi-waves model formulated in Ref. [23]. 

The raising, lowering, and population operators of the QD by denoted as a p q = 
!*%>,(/)(*%>,(?!> £md o Z p, q I ^p,g )(*%>,</1 \^p,g) {^p,g\i respectively. The total 
Hamiltonian is 


H = H d + H df + H t + AH, 
where the term 

fr _ ~ 

d ~ ^ a zp,q 

p,g 

describes the free electron motion, while the one 

Hj f = -(//E 0 ) ^ exp (*(p0! + # 2 )) + H.c. 
p,q 


( 3 ) 

( 4 ) 


( 5 ) 
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(H.c. stands for Hermitian-conjugate operator) accounts for the interaction of QD array 
with the electromagnetic field in the dipole approximation without tunnelling. Here, 
phases 


0i = (ka/2) cos (6), 02 = (ka/2) sin ( 0 ). 


( 6 ) 


represent the field delay per lattice period due to the finite propagation speed. The 
limit of 0i 2 —» 0 corresponds to the dense lattice, with ka y>> 1. 

The third term in Eq. (3) corresponds to the inter-dot coupling through the 
tunnelling. To interpret, it we introduce tunnelling coupling factors q ^ which are 
step functions of the discrete spatial coordinates, p and q: 


c( a , b ) 

^p,q 


£ (a ’ 6) , P, q = ~(N 1>2 /2 ),... - 1,0,1,..., N 1>2 /2, 
0, in all other cases. 


( 7 ) 


tion, 


The Hamiltonian of the tunnelling interaction is, in the tight-binding approxima¬ 
tor — — h ^ { a p+l,q\ + €pj\ a p,q) ( a p-l,q\ + ^p,q+l I a p,q) { a p,q+l I + £p,q\ a P,q) ( C 


p,q 


«EKSi.Jv«)< i >p+i,«i+Civ,><Vi,,i+C + iiW(Wii+Civ ? ><v ? -.i)' 


(b) 


p,q 


Coefficients in Eq. (8) are accounted for the tunnelling coupling between adjacent 
QDs in the excited and ground states, respectively. They are defined as scalars due 
to the isotropy of individual QDs and the array. The step-like structure in Eq. (7) 
corresponds to the operator notation for the tunnelling in the finite-size array: the 
number of the adjacent QDs coupled by the tunnelling reduces up to three for sites at 
edges of the array, and to two at corners. 

The last term in Hamiltonian (3) represents the local-field effects, in the Hartree- 
Fock-Bogoliubov approximation [17], [27]: 

A H = yN at gp a pg ( a p,q( a p,q) + K,q( a p,q)) ’ ( 9 ) 

p,q 


where p a jj and N a> p are components of the dipolar-moment vector and depolarization 
tensor of the single QD, V is the volume of the single QD, and angle brackets denote 
averaging of the corresponding operator with respect to the given quantum state. The 
summations over doubly repeated indices have been omitted in accordance with the 
usual convention. The depolarization tensor depends both on the QD configuration and 
quantum properties of electron-hole pairs, given by 

N a ,p = ^ J^J^{Y)\ 2 \^(Y l )\ 2 G a ^{Y-Y l )dfrdh l , ( 10 ) 

where £(r) is the envelope of the wave-function of electron-hole pair, and G a> p( r — r') 
is the Green’s tensor of the Maxwell’s equations in the quasi-static limit [14]. The 
present formulation may be applied to ID and 2D setting alike [24], Using the Hartree- 
Fock-Bogoliubov approximation, local-field effects for the single QD, were investigated 
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in the strong coupling regime in Refs. [17], and [27], and thereafter demonstrated 
experimentally [44]. 

The temporal evolution of single-particle excitations is governed by the Schrodinger 
equation, 

= #|#>. ( 11 ) 

The unknown wave function is taken in the form of a coherent superposition, 


|T) = V(^(f) e h/2)W i +# 2 -t)|op,,(t) e -h/2)Wi+# 2 -^)|^ ( 12 ) 


p,q 


where x l l pq (t), & p ,q(t) are probability amplitudes to be found. Projection of the 
Schrodinger equation onto the chosen basis leads to the following system of coupled 
nonlinear evolution equations for the probability amplitudes: 


d T 


M = 

dt p ’ q 




W^p,q Act) | &p t q ^tp,qi 




_M = _ 
dt p ’ q 


W* 1 + C«V ,-ie^ + Cn«W 


+ 1 

- ig^ P:q - «Aw|T P i9 | 2 $ m , 


— UP2 


(13) 


where F is the detuning parameter [27], g = — fiE 0 /(2h) is QD-field coupling factor, 
and Ac o = 4 tt/ ///,, 5 A r 0 : 5 . 3 /(/rV") is the depolarization shift. The normalization condition 
for system (13) is 

ATi/2 N 2 /2 

E E (| 1 - p ,,|= + |® p ,,p) = 1 . (14) 

p=—Ni/2 q=-N 2 /2 

The observable value of the energy is given by 

iVi/2 N 2 /2 

s = = - y. E [--f’d^i 2 - i^i 2 ) 


p=-7Vi/2 q=—N 2 /2 


— $ 


p,q \^p,q p 


- T 


p,q \^p 


n P i<pi _i_ c9>) <r> ,, P -i<pi _i_ /K 6 )<r> , p i<p 2 1 t( 6 ) <r> , 

-i,g e u- sp+p^p+i,*? 6 + u- sp } g+i^p,g+i e y 

a) T 1 \I/ ,1 e icpl ip "^ 2 4 - £ (a) J 


+9^pn^P,q + ^ |$p,g| + C.C. 


( 15 ) 


The analysis of the Rabi solitons in the 2D QD lattices is developed below on the basis 
of Eq. (13). 

In the simplest case, one may assume the inter-site coupling coefficients for the 
ground and excited states to be equal, = £ (h> = £. Next, we set, by means of an 
obvious rescaling, g = — 1 and sign (An.’) = — 1 . These signs imply the attractive onsite 
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linear interaction between fields and d> p (/ , and the self-focusing sign of the XPM 
(cross-phase-modulation, i.e., the cubic interaction between the different components) 
onsite nonlinearity. Actually, if g is originally positive, it can be made negative by 
substitution = — <f> M . And if Aui is originally positive, it may be made negative by 
means of the usual staggering substitution [6]. Thus, in Eqs. (13) there remain three 
independent parameters which can be combined into the frequency detuning, F, and 
the complex lattice coupling, £ exp (i0i j2 ). We here restrict to the basic interval, 
0 < 0i,2 A 7r/2, therefore both positive and negative values of F should be considered. 


2.2. The cylindrical-wave drive 


The spatial structure of Rabi waves is tunable by selecting the corresponding driving 
field. Of particular interest is cylindrical drive, defined as 

E(r, t) = e^Re |e 0 H^\k\r - r 0 |) exp (-hut) j , (16) 

where H^\x) is the Hankel function of the zero order and the first type. Such field 
may be excited by the infinitely thin and infinitely long current wire placed at point 
r 0 normally to the array’s plane (see Fig. 1 (b)). The current may be produced, for 
example, by a semiconductor quantum wire excited by the exciton-polariton [45]. The 
current point r 0 is placed at the center of the cell consisting of dots with the indices 
(p = 0,q = 0), (p = — 1,<7 = —1 ),(p = 0, q = —1) and (p = — 1,(7 = 0). The general 
statement of the problem is similar to the plane-wave case, the difference being only in 
the form of the H^f component of Hamiltonian (1), which is, now, 

Hdf h^z-E'O ^ ] ^p,q^P,q T H.C, (17) 

p,q 


where 


H, 


(i) 


S t 


+ 1/2) 2 + (g + l/2) 2 ) 


pq 


H S 11 (ka/V2 ) 


(18) 


The corresponding equations of motion are 
<9T r 




dt 

+i 


= iF^ 


p,q 


£( a )\I/ -i 4 - F a \ C P+1,9 4 - £( a )\Er lC M-i I tW 
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( 20 ) 
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Ni/2 JV 2 /2 

E E (i*pJ 2 +l* M i 2 ) is m i = i- (21) 

p=—Ni/2 q=-N 2 /2 

Actually, Eq. (19) describes other types of driving field too. Namely, substitution 
k —>■ ik (for real k) transforms the Hankel function into the correspondent Macdonald 
function Kq(x), which represents driving the Rabi wave by surface plasmons guided by 
a noble-metal wire [14], or carbon nanotube [46]. 


3. 2D Rabi solutions 

3.1. Dispersion relations 

In the limit of A u> = 0, Eq. (13) without the local-field terms reduces to the 2D 
generalization of the Rabi-wave model considered in Refs. [ 22 , 24], In this case, the 
dispersion relation can be derived by looking for solutions as 

{^ P , q , = {A,B}exp(i(K x p + Kyq))exp(-i£lt), (22) 

where n x>y are wavenumbers, A, B are unknown wave amplitudes, and Q is an unknown 
frequency. 

A straightforward analysis yields two branches of the dispersion relation, 

Di )2 (kz, K y ) — - 2£ [cos (k x ) cos (0i) + cos (k v ) cos (0 2 )] 

± yjg 2 + {F - 2£ [sin (k x ) sin (0i) + sin (Ky) sin (0 2 )]} 2 , (23) 

examples of which are shown in Fig. 2 for g = — 1 and £ = 1. 

In the limit of <pi = <f 2 = 4> = 0, Eq. (23) goes over into the known dispersion 
relation for the usual system of linearly coupled 2 D DNLS equations, cf. Ref. [ 6 ]. It 
features two similar branches shifted by a constant, A Q = 2 g 2 + F 2 , and corresponds 
to the Rabi waves. In the case of cf) 1 , 02 0, shift AD is no longer constant, as it 

depends on wavenumbers k x and K y . 

Note that 

fi(«*, Ky, 01, 02) = fi(-K®, -Ky, -01, -02), (24) 

but 

fl(, Ky, 0i, 02) 7 ^ K x , K y , 0i, 0 2 )- (25) 

Equation (24) shows that the inversion of the propagation direction of the driving 
field inverts the direction of the Rabi-wave propagation. This is according to the non¬ 
reciprocity of Rabi-waves in ID chains, noted in Ref. [39]. 
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la- 








Figure 2. (Color online) Dispersion curves $7 = 0,2 (/%,, k 9 ) for the linearized system 
with the plane-wave excitation: F = 0 (black surfaces), F = 1 (red surfaces), and 
<j>i = (j >2 = 0 (a), <j>i = (j >2 = 7t/4 (b), (j >i = 0, (j >2 = y/Zit/A. (c). Other coefficients are 
£ = 1 and g = — 1. Discrete solitons are expected to exist in the semi-infinite gap. 


3.2. Constructing Rabi solitons and their radiation patterns 


Equations (13) may give rise to stationary discrete 2D fundamental solitons and vortices, 
looking for localized solutions as 


vh = 

* p,q 


—iQt 


A 


-pm 


$ = 


—iQt 


B 


p,qi 


(26) 


where 0 is the corresponding carrier frequency, and A Ptq , B p q are localized complex 
lattice helds vanishing at infinity (in terms of the numerical solutions, they vanish at 
boundaries of the computation domain). Stationary soliton solutions of Eqs. (13) were 
numerically obtained by adopting the nonlinear equation solver based on the Powell 
method while direct dynamical simulations were based on the Runge-Kutta procedure 
of the sixth order, [47], [48]. The numerical solution was constructed for a finite lattice 
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of size 21 x 21 (unless stated otherwise), with periodic boundary conditions. We have 
checked that a larger lattice produces the same results. 

The emitted held is described by the held operator in the Heisenberg representation 

[ 15 ], 

i_i d\ r i 


E(r ' t) = -4^ (V ' V -^ 


lv 


r — r 


7^~ 


r — r 


)dV. 


(27) 


Here, operator P of the induced polarization describes the displacement current induced 
by the soliton’s prohle, which may be considered as an external source, from point of 
view of the antenna. The onsite polarization operator, written in terms of creation- 
annihilation operators is 

P P ,q(t) = pa+ q (t) + H.c. (28) 

Due to the smallness of the held produced by the soliton in the QD array, its emission is 
determined by the linear approximation, hence the total held is built as a superposition 
of partial helds emitted by different QDs independently. The plane-wave drive is, thus, 

V 


E(r ,t)= - 


47Te 0 


V ' V -?Jp) e ' 


X 


M/2 


M/2 

^ E 

p=-M/2 <Z=—m/2 


Pzp,q (p c |r e x pa Gyqof) 
|r - e x pa - e y qa\ 


(29) 


To obtain the observable value of held, it is necessary to replace the polarization operator 
by its main value with respect to the quantum state (12), (T|P 2 Pi9 |T). The result is 

P zp ,q(t ) = ( P sp ,q(t )) = ^T P;? (f)<f>* g exp [i{pct>i + qfa)] exp (-iut)+ c.c.(30) 

It is convenient to present the far-zone held in the spherical coordinates with the origin 
set at the central point of the cell with p = q = 0, dehned as x = R sin ( 9 ) cos (p), y = 
R sin ( 0 ) sin (p), z = R cos ( 6 ). The terms of order 0(R~ 2 ) should be omitted, and the 
radial (longitudinal) component of the electric held vanishes too. As a result, we obtain, 
for the quasi-spherical observable held: 

iT 


E 


Rad 


= lim E 


pe e of 1 e-^h- R / c ) 


R —^oo 

with the radiation pattern 


47reo c 2 


R N 


R 


F e,<p,<f>i,<f>2;t- 


T c.c., (31) 


F(6,p,f> 1 ,f> 2 ]t -) = sin# 

c 

f \UJCL . - 

x exp Up — sm 6 cos ip + 20i 
L c 


Ni/2 N 2 /2 

Y, E 


p=-N i/2 q=-N 2 /2 
uj a 


+ iq 


c 


sin#sin p + 20 2 }, 


(32) 


where t = t — R/c + c~ l sin # \p cos p + q sin p]. We stress that, in contrast with classical 
macroscopic antennas, the radiation pattern given by Eq. (32) is non-steady, as it 
exhibits a slow dependence on the spatiotemporal variable, t — R/c, which corresponds 
to the amplitude and frequency modulation of the held due to the time dependence of 












Soliton nanoantennas in two-dimensional arrays of quantum dots 


12 


probability amplitudes in quantum state ( 12 ). The stable soliton produces, through Eq. 
(32), a radiation pattern in the form of 

JVi/2 JV 2 /2 

F(9,p, 0i,0 2 ) — sin# ^2 \,qBp,q ex P Sin6*cos (p+20! j 


p=-N 1 /2q=-N 2 /2 

x exp \ iq 


sin 9 sin p + 20 2 J | • 


(33) 


In particular, the radiation pattern produced by the ID in the 2D space can be written 
in a simple form by setting IV 2 = 0 in Eq. (33): 

N/2 

F(9 , p, 0) = sin (9) A P R* ex P {W —— (sin #) ( cos a) +20j |, (34) 


p=—N/2 


where the notation is changed as A p/j —> A p , B Ptq —> B p , 0i -7 0, N± -7 N, to comply 
with that adopted in Ref. [27]. 

It is relevant to note that, according to Ref. [39] the radiation pattern is subject 
to the Onsager kinetic relations, which take into regard symmetry constraints [40]. As 
it follows from Eqs. (13) and (26), (A_ p _ q Bf p _ q ) = (A Ptq B* q ), hence 


F(tt - 9, p, -0 1 , -02) = F(0, p, 01, 02), 


(35) 


while F(tt — 9, p, 0i, 0 2 ) 7 ^ F(9, p , 0i, 0 2 ). It shows that, rotating the QD array by 180°, 
one should invert the propagation direction of the driving held to keep the physical state 
invariant, similar to the the sign of the angular velocity in the rotating liquid, or the 
sign of the magnetic held [40]. This symmetry agrees with the non-reciprocity of the 
Rabi waves exhibited by Eq. (24). 

For the cylindrical drive, the radiation pattern in the far-held zone is obtained from 
Eq. (31) in a similar way: 

JVi/2 N 2 /2 , . . 

F(9, p) = sin 9 E E Ap, q B* tq I S P) q I exp l -sin 9 [p cos p + q sin p\ 1 . (36) 

p=-N 1 /2q=-N 2 /2 ^ ' 

The useful information in antenna context are the radiation patterns in the E- and 
H -planes dehned, respectively, as F e {9) = F(9,p = n/2) and F/ l (6 l ) = F(9 = n/2,p) 

[28]. Plots of these quantities follow the corresponding dynamically stable soliton 
evolution in hgures below. 


4. Results and discussion 

4.1. Rabi solitons driven by the plane wave 

After stationary solitons were found as outlined above [see Eq. (26)], their stability 
against small perturbations was examined in the framework of linearized equations for 
small perturbations added to the stationary solution, 

= [Ap, q + SAp^it)] e lClt , $ Pl q(t) = [B Ptq + 8B p%q {t)\ e lClt . (37) 
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In this case, the linearization of Eq. (13) with respect to the perturbations SA P:q and 
SBpq, leads to a linear eigenvalue (EV) problem for the instability growth rate of small 
perturbations. Obviously, the soliton is stable if all the eigenvalues have zero real parts. 
On the other hand, the presence of EVs with a positive real part indicates at exponential 
or oscillatory instability, in the cases of purely real or complex EVs, respectively. Finally, 
the dynamical stability of the discrete solitons was verified by dint of direct simulations 
of Eqs. (13). 

It is well known that the single-component 2 D DNLS equation gives rise to three 
types of fundamental discrete solitons: stable onsite, and unstable hybrid and inter-site 
ones [49], along with several types of vortices with vorticity S = 1 (eight configurations) 
and S = 2 (four configurations) [49, 50, 51]. The vorticity is identified as the winding 
number of the phase of the lattice field, i.e., the total change of the phase along a closed 
curve surrounding the pivotal point of the vortex, divided by 2n. We here consider the 
values of S' = 1 and S = 2 . 

Obviously, the soliton complexes must be located inside the semi-infinite gap of 
dispersion relation (23), see Fig. 2. The most straightforward counterpart of the single¬ 
component 2 D DNLS equation is the model of the QD array based on Eq. (13) with 
zero detuning, F = 0 and ch = 62 = 0. In this system, we find soliton complexes 
formed by two identical real components of the onsite, hybrid, and inter-site types (Fig. 
3), and families of discrete vortex solitons with topological charges S = 1 and S — 2 
(together with quadrupoles), which are supported by the attractive XPM interaction. 
The localized modes in each component are similar to the corresponding types of 
fundamental and vortex solitons formed in the single-component 2 D DNLS equation. 

Different types of vortices in each component are schematically shown in tables I 
and II. Configurations a, c, d, e, f, g, i, and j represent possible structures of onsite 
vortices, with the pivotal point located on a lattice site, while b , h, k, and l generate 
off-site vortices, whose pivot falls located between lattice sites (the same nomenclature 
as in Ref. [50]). An example of the vortex of type a is plotted in Fig. 4. 

The linear-stability analysis shows that the complexes formed by identical 
fundamental-soliton onsite components are stable, and that there are stability regions 
in the respective parameter space for vortex configurations of types a, b (ground-state 
vortices), e, g with S — 1 (excited state vortices), as well as for types j, k with S = 2 . 
The corresponding radiation-field distributions for the onsite fundamental solitons are 
plotted in Fig. 5, while the same for the stable vortices of the a and b types is shown 
in Fig. 4. The vortices of type k actually feature an offsite quadrupole configuration 
(rather than a truly vortical one), based on the four-site frame and sharp phase changes 
in steps of 7 r between lattice sites. The stability regions of all the above-mentioned 
complexes almost overlap with the stability regions of their counterparts in the single¬ 
component 2 D DNLS equation. Briefly speaking, the vortex-stability range is narrow, 
particularly for spatially broad vortices, and shrinks with the increase of £ (which is the 
counterpart of the inter-site coupling constant C in single-component DNLS equation 
[50]). This is in accordance with the fact that large £ corresponds to the continuum 
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limit, in which all the localized vortex modes are unstable. We have also checked a 
possibility to form vortex complexes formed of two identical vortices with opposite signs 
(phase shift 7r), but they all turn out to be unstable. 

The predicted dynamical stability and instability has been confirmed by direct 
simulations for the onsite fundamental soliton complexes and certain vortex complexes 
for typical values of the system’s parameters. The onsite fundamental complexes and 
those vortex complexes which were predicted to be stable indeed keep their amplitude 
and phase structure (i.e., the vorticity) in the course of the evolution. On the other 
hand, hybrid and offsite complexes, whose instability is predicted by the linear-stability 
analysis, radiate away a significant part of their norm and rearrange into dynamically 
robust onsite-centered localized breathers. Direct simulations of vortex complexes, 
which were predicted to be unstable, demonstrate that their amplitude profiles are 
quite robust, and they keep the vorticity in the course of the evolution. However, 
the symmetry with respect to the pivot is broken, and the symmetry between their two 
components is broken too. The so established mode may be categorized as an irregularly 
oscillating vortex breather, see Fig. 6. 





Figure 3. Families of fundamental solitons, formed by two identical onsite (solid 
black line), hybrid (solid red line) and offsite (dotted green line) components, are 
represented by the corresponding |Aw(f2)| dependencies in (a). Examples of these 
solitons for $7 = —7 [the vertical blue line in (a) which intersects all three curves] are 
shown in (b) - onsite, (c) - hybrid and (d) - offsite, respectively. This figure pertains 
to the model with the plane-wave excitation, see Eq. (13). 

In the 2D array of QDs with nonzero detuning (F = 1 is chosen below, as a 
characteristic example), all types of the above-mentioned fundamental and vortical 
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type a, S=1 



o 



type b, S=1 





Figure 4. Amplitude, phase profiles, and radiative patterns in E- and ff-planes of 
vortices with 5 = 1 defined as types a and b in Ref. [51]. The parameters in Eq. (13) 
are £ = 0.01, F = 0, <j>\ = = 0, andfl = —19.3. Both are stable according to the 
linear-stability analysis and dynamical simulations. 


complexes can be generated, but the symmetry between the two components {A pq and 
Bp q} is not preserved. Simultaneously, the similarity to the single-component DNLS 
equation is lost. Namely, the overall shape of the components in fields A pq and B pq is 
kept, but amplitudes and widths of the components change, and their precise symmetry 
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Figure 5. The radiative patterns in the E- and iJ-directions, F e (9) and Fh(4>), for 
the complex consisting of two fundamental onsite solitons with (a) fl = —7 , F = 
0, 4>\ = </>2 = 0 and (b) $7 = —7, F = 0, 4>\ = </>2 = 7t/ 4. The plane-wave excitation is 
considered, see Eq. (13). 


about the central point is broken (slightly violated). In this case, the linear-stability 
analysis does not predict the existence of genuinely stable fundamental complexes, with 
the exception of the vortex complexes of type a, with S — 1 , and of type k, with S = 2, 
see below. However, dynamical simulations demonstrate that localized modes of the 
onsite type form robust breathing complexes, which keep the initial overall structure (in 
particular, the vorticity is kept). 

On the other hand, the computation of the eigenvalues for small perturbations 
shows stability of vortex complexes of type a, with S — 1 , as well as of types k with 
S = 2 in the presence of F 0. Dynamical simulations confirm their stability in the 
corresponding parameter regimes. Other families of vortex complexes formed at F 0 
are unstable according to both the linear-stability analysis and dynamical simulations. 
In fact, the instability destroys their phase structure, while the amplitude patterns 
remain quite robust. 

Thus, the general conclusion concerning the system (13) driven by the plane wave is 
that, in the absence of the phase shifts, <f\ = (p 2 = 6, the the stability of the complexes is 
similar to what is known about their counterparts in the single 2D DNLS equation single 
2 D lattice case. Namely, irrespective of the presence of the detuning, F, in Eq. (13), 
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Figure 6. The evolution of one component of an unstable vortex of type c with 5 = 1. 
(a) The amplitude and phase profiles of the stationary vortex with fl = —10. Panels (b) 
and (c) show the profiles in the course of the perturbed evolution (t = 50 and t = 100 
in arbitrary units, respectively). Other parameters are £ = 0.01, F = 0, d>i = 02 = 0. 


the those modes which were stable in the single-component model, remain completely or 
effectively stable in the two-component one, evolving in some cases from the stationary 
shapes into breathing ones, but keeping the overall structure. 

Considering the system with non-zero phase shifts, 0i and 02 in Eq. (13), we mainly 
focus on two configurations, viz ., the diagonal and anisotropic ones, with F\ = O 2 = tt/4 
or = 0, 02 = \/27r/4, respectively. In both these cases, we have found complexes of 
all the above-mentioned types. Namely, fundamental solitons of the onsite, hybrid and 
inter-site types, with symmetric components, in the case of F = 0 (see Fig. 7), and 
their asymmetric counterparts at F ^ 0. The existence boundary for the fundamental 
soliton complexes is slightly altered by the nonvanishing phase shifts, as can be seen 
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from small changes of the respective dispersion surfaces displayed in Fig. 2 . Example 
of the soliton radiation fields for a stable configuration is shown in Fig. 5(b). 

Counterparts of all the vortex field configurations, which are identified above in 
the system with </>i j2 = 0, are also found in the system with bi .2 0. Differences are 
observed in the boundary of the existence region, and parameter regions where robust 
vortex breathers appear. In Fig. 8 , the evolution of the symmetric vortex soliton 
of type a is displayed for <fi = </> 2 = 7 t/4 , F = 0. According to the computation 
of perturbation eigenvalues, stable complexes of two a type vortices (with topological 
charge S = 1 ) keep their stability at hi = <b 2 7 ^ 0 . while unequal phase shifts destabilize 
the complexes. The radiation pattern keeps its shape, as shown in Fig. 8 (c). Most 
robust are complexes formed by the quadrupoles of type k with S = 2 , in the sense 
that they are stable in a certain parameter area for all relations between phases Si and 
ct> 2 • Lastly, the phase differences always destabilize other types of vortex complexes. In 
general, the value of the corresponding instability rate, which is proportional to a purely 
real eigenvalue, increases with increasing <pi, 02 and fixed F. Direct simulations confirm 
that only complexes formed by the onsite fundamental solitons, and certain types of 
vortices are stable in particular parameter ranges, chiefly in the form of robust localized 
breathing patterns. Similar to what is mentioned above, the instability destroys the 
phase structure of vortices, without essentially affecting the amplitude pattern. 

Thus, we conclude that the changes in values of phase shifts Si )2 i n Eq. (13), which 
are introduced by the oblique incident plane waves, as well as the detuning, F, do not 
cause qualitative changes in the structure and effective stability of the 2D self-trapped 
fundamental modes and particular types of vortex complexes, of type a with S = 1 , and 
type k with S = 2 (the quadrupole). Finally, systematic simulations demonstrate that 
moving discrete-soliton complexes cannot be found in the present 2D model. 

4-2. Rabi solitons driven by the cylindrical wave 

For the 2D QD array with the cylindrical-wave excitation, described by Eq. (19), we 
report here results obtained for ka = 7 t/ 4 . As we shown by detailed analysis, these 
results adequately represent the generic situation. 

Similar to the planar-wave model, soliton complexes are located in the semi-infinite 
gap in the linear spectrum. However, unlike Eq. (23), the spectrum with the cylindrical 
driving wave cannot be found analytically, therefore it was calculated numerically. Two- 
component fundamental solitons are built, as above, of two components, each of them 
being of the onsite, inter-site, or hybrid type. Vortex complexes found in the present 
system are displayed in Fig. 9. These are complexes of b and l types with S' = 1 , or of 
k and h types with S = 2 , see Tables I and II. In the system with a nonzero detuning 
parameter F^O, all types of the above-mentioned fundamental and vortical complexes 
are generated, but the symmetry between the two components is not preserved. 

The linear-stability analysis shows that the complexes formed identical 
fundamental-soliton onsite components are stable for all value of detuning F, and there 
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Table 1 . Schemes of the amplitude and phase profiles corresponding to the two 
identical components of the discrete vortex with S' = 1 in the system with F = 0, 0i = 
02 = 0 [the model with the plane-wave excitation, Eq. (13)]. Explicitly written shares 
of 7r (or 0) are values of the phase at main sites carrying the vortex complex. Symbol 
x designates sites with zero amplitude. Only complexes with the components of the b 
and l types for S = 1, and of the k and h types for S = 2, respectively, are found in 
the model with the cylindrical-wave excitation. 


Type 

S 

Phase configuration 

Stable (plane wave/ cylindrical wave) 

a 

1 

0 

7t/ 2 x 37 t/2 

7T 

Yes/No 

b 

1 

0 3tt/2 

7t/2 7T 

Yes/No 

c 

1 

7t/4 0 77t/4 

7t/ 2 x 37 t/2 

37t/ 4 7r 57 t/4 

No/No 

d 

1 

0 

7t/ 4 x 77 t/4 

7t/ 2 xxx 37 t/2 
37t/ 4 x 57 t/4 

7r 

No/No 

e 

1 

0 

7t/ 4 x 77 t/4 

7t/ 2 x 0 x 37 t/2 

37t/ 4 x 57 t/4 

7T 

Yes/No 

f 

1 

0 

7t/4 0 77t/4 

7t/2 7t/ 2 x 37 t/2 37t/2 

37t/4 7t 57t/4 

7T 

No/No 

9 

1 

0 

7t/4 7T 77t/4 

7t/2 37t/ 2 x 7r/2 37 t/2 
37t/ 4 0 57r/4 

7T 

Yes/No 

h 

1 

0 77t/4 

7t/4 0 37t/2 37t/2 

7t/2 7t/2 7T 57t/4 

37t/4 7T 

No/No 
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Figure 7. Examples of fundamental-soliton complexes [in the model with the plane- 
wave excitation, Eq. (13)], formed by identical components, at F = 0, ft = —7, and 
fi = 02 = tt/ 4. Plotted are real and imaginary parts of the stationary lattice fields 
for the complexes of the onsite (a), hybrid (b), and intersite (c) types, respectively. 


are stability regions in the respective parameter space for vortex configuration of type 
k with S = 2 at F = 0. The stability region of latter mode almost overlaps with its 
counterpart for the corresponding vortex complex in the 2D system with the plane- 
wave excitation. Radiation patterns of stable solitary Rabi structures excited by the 
cylindrical wave are plotted in Fig. 11. 

The predicted dynamical stability has been confirmed by direct simulations for 
the onsite fundamental soliton complexes (Fig. 10) and vortex complexes for particular 
values of the system parameters. Similar to what was seen above, stable vortices preserve 
their amplitude and phase patterns, as shown in Fig. 12. In general, effectively stable 
localized patterns (including those which are unstable in terms of the eigenvalues but 
exhibit persistence in direct simulations) evolve as breathing modes. In contrast, the 
unstable hybrid and inter-site complexes radiate into background a significant part 
of their energy. The remaining energy can be reorganized as a new localized onsite 
breathing structure, see Fig. 13. On the other hand, unstable vortex complexes are 
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o 




p 




Figure 8. Initial (a) amplitude and phase profiles of the vortex of type a, in terms of 
Ref. [51], the result of its evolution at t = 100 (b), and the respective radiative patterns 
F e and (c). The parameters in Eq. (13) are £ = 0.01, F = 0, = fa = 7t/ 4, and 

n = -is. 


more robust structures with significantly reduced radiation of energy to the background. 
The newly formed breathing structures preserve localization but phase coherence is 
destroyed. The symmetry of the profile with respect to the center of vortex is broken, 
as well as the symmetry between the components, see Fig. 14. 

In the physical system underlying the present lattice model, parameters may be 
known only with a finite accuracy, which makes it necessary to test the stability of the 
localized modes with respect to random variations of coefficients of Eqs. (13) and (19). 
We have performed this test by introducing random changes, of relative amplitude 
10%, of coefficients £ and F, which represent uncertainties in QD energy levels, and 
tunneling matrix elements that determine the strength of the inter-site coupling. The 
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Table 2. Schemes of the amplitude and phase profiles corresponding to vortex solitons 


with 5 = 2. 


Type 

s 

Phase configuration 

Stable (plane wave/cylindrical wave) 

i 

2 

tt/2 0 3tt/2 

7T X 7T 

3tt/2 0 tt/2 

No/No 

j 

2 

\ x \ 

to to 

o X X X o 

CO 

to to 

Yes/No 

k 

2 

0 7T 

7T 0 

Yes/Yes 

i 

2 

0 3vr/2 

7t/ 2 0 7T 7r 

7T 7T 0 7t/2 

3tt/2 0 

No/No 




Figure 9. The amplitude (a) and phase (b) patterns of the vortex in each component 
of the complex of the k type with S = 2, excited by cylindrical wave. 


result (not shown here in detail) is that all the localized modes which were found above 
to be dynamically stable are also structurally stable against the random variations of 
the system’s parameters. In fact, the robustness of discrete solitons against random 
deformations of the underlying lattices is known in earlier studied models [6]. 

The bottom line of this subsection is that in the system with the cylindrical-wave 
excitation, fundamental onsite soliton complexes and vortical ones of the k type with 
S = 2 are stable in certain areas of their existence region. In addition, we point out that 
the soliton-produced radiation patterns may be used as digital signals for data coding 
in the framework of the quantum information transition. Due to the different patterns 
of these signals for different types of the models, driven by the plane and cylindrical 
waves, they may be identified clearly. 
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Figure 10. The evolution of the amplitude pattern of the onsite soliton complex with 
// = —17, A to = —30, F = 1 in the model with the cylindrical-wave excitation, see Eq. 
(19). Plots correspond to t = 0, 50, 100 presented in arbitrary units. 


4-3. The concept of solitonic nano-antennas 

As said above, the transmitting antenna is defined as a device which converts the 
near field into the far field [28]. On the other hand, the Rabi solitons feature the 
spatially distributed polarization induced by the RO, as seen in Eqs. (28) and (30). 
The 2D electromagnetic radiation patterns induced by these polarization profile suggest 
to consider the present system in the context of the realization of transmitting antennas. 
Thus, the Rabi soliton offers a previously unexplored physical mechanism to emulate the 
antennas, while the driving field plays the role of the external energy source launching 
and driving this device. 

For the data transfer, the emitted field should be temporally modulated by an input 
signal through the driving electromagnetic field. Actually, this modulation is relatively 
slow, hence it may be considered as an adiabatic process. As mentioned above, a basic 
characteristic of the antenna is the 2D radiation pattern, or, in some cases, partial 
patterns in the E- and H- planes [28]. The necessity to consider the radiation profiles 
in 2D makes the present analysis principally different from the study of RO solitons in 
ID arrays, which was recently reported in Ref. [27]. 

Radiation pattern (33) produced by the QD array can be rewritten in terms of the 
standard notation adopted in the antenna theory [28] as F(6, <p) = sin ( 0 ) • (AF)n iX n 2 . 
The first factor, sin(0), is the radiation pattern of a single emitter (in our case, the 
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Figure 11. Radiation patterns F e and Fh for (a) the onsite complex with // = 
— 17, A to = —30, F = 1, and (b) the vortical complex of the k type with F = 0,p = 
17.2, S = 2, in the model with the cylindrical-wave excitation, Eq. (19). 


electric dipole), while the second one is the array factor, 

N\/2 N 2 /2 

(AF) NixN2 = y y: A p,q B p, q exp {ip [c^uja (sin 9) (cos <p) +20j } 

p=—N\/2 q=-N 2 /2 

x exp {zg [c _1 cna (sind) (sincp) +2</> 2 ] } . (38) 

In the macroscopic phased antenna array, elements are excited by a signal having a 
constant amplitude and the phase which uniformly increases along the array’s axis [28]. 
The array factor for that case corresponds to Eq. (38) with A p ^ q B* q = 1. The phase 
variation allows one to perform scanning of the space by turning off the main lobe of the 
radiation pattern. One of the most efficient means to implement the antenna control is 
the use of switched-beam systems which can impose different angular patterns, in order 
to enhance the emitted signal in a preferable direction. The beam-forming algorithm is 
implemented through a complex pattern of the excitation of individual elements [which 
corresponds to factor S pq = A P:q B pq in notation (38)], adjusted so as to place the 
maximum of the main beam in the necessary direction [28] . 

The purport of the solitonic nano-antenna concept is to introduce the beam-forming 
algorithm, which is determined by the set of factors S pq , using the soliton-emission 
profile. As several different stable solitons may exist in the given array, the choice of 
the overall profile depends on initial conditions. This implies that a finite number of 
predefined array factors exist for a given set of antenna parameters, while the initial 
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Figure 12. The evolution of the stable vortex complex of the k type, excited by 
the cylindrical wave, is illustrated by amplitude and phase patterns, shown at two 
moments of time. Parameters are F = 0, y = 17.2, S' = 2. Plots correspond to t = 0 
and 50 in arbitrary units. 


conditions give a facility for choosing and switching the suitable one. To the best of our 
knowledge, this feature has no analogs in previously developed types of nano-antennas, 
and it seems quite promising from the practical point of view. 

For the stable solitons built of two identical components, A and B in terms of Eq. 
(22), the array-factor pattern is independent on the phase profile of the soliton. Partial 
radiation patterns in the E- and //-planes for different types of the solitons and different 
types of the excitation are displayed in Figs. 4, 5,8, and 11. 

A basic problem that one should take care of to implement emitting nano-arrays 
in the real physical setting is maintaining the coherence of individual emitters in the 
array over sufficiently long time. The destruction of the coherence is contributed to by 
physical mechanisms of dephasing and relaxation, as well as by material imperfections 
of the underlying array. Recent experimental results clearly demonstrate that these 
problems may be resolved, allowing quite large QD arrays to maintain the coherent 
behavior even at room temperatures [52, 53]. 

There are two types of QDs that are commonly employed in nanophotonics [54] . The 
first one is built of semiconductor nanocrystals (typically II-IV compounds) embedded 
in a glass matrix. In Ref. [52], the QD ensemble of the CdSe/ZnS core/shell 
nanocrystals, capped with octade-cylamine, has been produced. It was fabricated as a 
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Figure 13. The perturbed evolution of the amplitude pattern of an offsite soliton, for 
jx = —17, Aw = —110, F = 1. The cylindrical-wave excitation is used, see Eqw. (19). 
Plots correspond to t = 0, 50, 100 in arbitrary units. 


2D inhomogeneously broadened close-packed network of 50 x 50 sites, with the average 
diameter and interdot distance 5.2 and 7.9nm, respectively. There also appeared dark 
QDs associated with defects in the ensemble. The qualitative picture of the exciton 
dynamics in the QD array observed in Ref. [52], based on direct measurements, shows 
that when an exciton is photoexcited in a high-energy QD, energy transfer occurs 
preferentially to a low-energy QD. At low temperatures, the exciton is trapped at a 
local low-energy site to which the energy was transferred. In contrast to that, at room 
temperature the exciton hops repeatedly, until it is transferred to a dark QD, where 
it undergoes nonradiative recombination. When an exciton is photoexcited in a low- 
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Figure 14. The evolution of an unstable vortex of the l type in the model with the 
cylindrical-wave excitation. Parameters are F = 0,y = 17.2, S' = 1. Plots correspond 
to t = 0 and t = 100 in arbitrary units. 


energy QD, it tends to be trapped in it. However, at room temperature, there is a 
non-negligible probability for the exciton to be transferred from a lower-energy QD to 
a higher-energy one. In Ref. [52], time- and spectrally resolved fluorescence intensities 
were measured by means of site-selective spectroscopy at both room and low (80 K) 
temperatures. The respective inverse radiative decay rate is found to be 15 ns, which 
makes such type of arrays promising for the implementation of nano-antennas. 

The second type of QDs is the self-organized structures produced by the epitaxial 
crystal growth in the Stranski-Krastanov regime [54], The self-assembled QD array was 
grown in Ref. [53] with the help of the molecular-beam epitaxy on a semiconducting 
(lOO)-oriented GaAs background, with a 500 nm N + -GaAs buffer layer. The respective 
lens-shaped Ino. 65 Alo. 35 As quantum dot is a part of a sphere with a fixed height of 3.4 
nm and base diameter of 38 nm. The photoluminescence spectrum of the QD array 
is time-resolved at the excitation density of 1065 W/cm 2 , at 77 K. The corresponding 
exciton lifetime, found from the measurements, is 800 — 1200 ps, which complies with 
the respective value for the arrays of the first type [52], Note that, for the second type, 
the growth of the top layer allows the cavities to be formed and the electrical contacts to 
be applied, which makes it especially suitable for the implementation of nano-antennas. 

Lastly, it is relevant to mention that the soliton mechanism of the self-trapping 
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of the antenna helps to mitigate detrimental effects of imperfections of the underlying 
QD lattice, as discrete solitons are well known to be robust against deformations of 
the lattice [6]. Moreover, the fact that the phase structure of topologically organized 
solitons, such as discrete vortices [51], is stable in imperfect lattices, also indicates that 
the so constructed nano-antennas tend to stabilize themselves against dephasing. 

It is relevant to note that the self-organized lattice built of semiconductor QDs not 
a unique structure allowing the implementation of soliton-based nano-antennas, and the 
RO is not a unique enabling mechanism for this. Another promising way for achieving 
this purpose is suggested by theoretical analysis of discrete dissipative plasmon solitons 
in an array of graphene QDs [55]. The single QD in such an array represents a doped 
graphene nanodisk placed on top of the plane background. The single-dot excitation 
represents a confined surface plasmon with the resonant frequency in the THz or infrared 
range. The QDs in the array are coupled by long-range dipole-dipole interactions. As 
it was demonstrated in [55], the soliton formation and its stability take place under the 
control of an incident driving electromagnetic wave via the Kerr optical nonlinearity. 
In spite of the different physical origins, both the surface plasmon solitons predicted 
in [55] and Rabi solitons considered above exhibit strongly confined one- or two-peak 
areas of electrical polarization. Such peaks may be symmetric or asymmetric, while the 
spatial structure of the confinement area is tunable via the incidence angle of the oblique 
driving field. Thus, the above-mentioned graphene QD-array makes it possible to design 
an electrically controlled nano-antenna for THz and infrared frequency ranges. The 
radiation pattern (or array factor) for this antenna is given by Eq. (38) with necessary 
modifications (the slowly varying amplitudes of orthogonally directed dipole moments 
can be found as a solutions of the nonlinear coupled equations derived in [55]). The 
qualitative shape of such a radiation pattern is similar to that of the Rabi-soliton-based 
nano-antenna proposed above. 

5. Conclusions 

The objective of this work is to introduce the concept of the tunable nano-antenna 
array based on the discrete-soliton patterns formed in the 2D nonlinear lattices of 
semiconductor QDs (quantum dots). These lattices can be realized as square-shaped 
arrays of identical two-level quantum oscillators (the self-organized semiconductor QDs), 
coupled to nearest neighbors by the electron-hope tunneling and interacting with 
the external electromagnetic field. The local-field corrections, which account for the 
difference between the field inside the QD and the external field, induce the nonlinearity 
of the electron-hole motion inside each QD. 

The main conclusions of our study are summarized as follows. 

(i) The model of RO (Rabi oscillations) in the 2D QD lattice has been derived, 
taking the inter-dot tunnelling and local-field correction into account. The model is 
based on a set of linearly and nonlinearity coupled DNLS equations for probability 
amplitudes of the ground and first excited states of two-level oscillators (QDs). 
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Two different driving electromagnetic fields were considered, viz ., the plane-wave and 
cylindrical ones. In the former case, the coupling coefficients are complex, with absolute 
values independent of the QD position, and phases linearly increasing across the QD 
array in both directions. For the cylindrical-wave drive, the absolute values of the 
coupling coefficients depend on the distance between the given QD and the source of 
the driving field. The corresponding phases depend nonlinearly on the distance. 

(ii) Stable discrete-soliton complexes are found. They are built of onsite 
fundamental single-peaked solitons, taken in both components, or discrete vortex 
solitons of certain types. 

(iii) The emission properties of stable solitary modes have been characterized by 
angular radiation patterns. These patterns strongly depend on the type of the discrete 
localized mode. 

(iv) The concept of self-assembling nano-antennas, based on the stable discrete- 
soliton complexes in the nonlinear lattices, is introduced. The necessary type of the 
localized mode may be selected by the initial conditions, which, in turn, can be controlled 
by the external optical field (to be supplied in the form of a strong laser pulse). As a 
consequence, a finite number of predefined radiation patterns can be provided by a 
given antenna. This way of the operational control of nano-antennas has no analogs in 
previously developed antenna schemes. The stability of the self-trapped nano-antennas 
against structural imperfections and intrinsic dephasing has been considered. Thus, the 
system proposed here can be related to switched-beam systems used in macroscopic 
antennas [28], which is a promising setting for applications to nanoelectronics and 
nanooptics. 

(v) The concept of soliton-excited nano-antennas may be carried over to the surface- 
plasmon mechanism of the soliton formation in the array of graphene QDs with the 
Kerr nonlinearity predicted in [55]. This type of the nano-antennas is promising for 
applications in the THz and infrared frequency ranges. 
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